A spatial network explanation for a hierarchy of urban power laws 
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In this Letter we report an application of spatial complex networks to simulate the evolution 
of the spatial economic system. The model depicts the spatial trade network with land areas as 
nodes and trade relations between activities as connections. The network grows by addition of 
new pairs of activities using simple market economic assumptions. The definition of nodes and 
connections allows us to compare model and reality on the level of market land values. The model 
provides an explanation to several empirically observed statistical phenomena in spatial economics. 
Furthermore, market demand is an important driving force for most other spatially extended systems 
of human origin. Thus, a model that simulates the geographic distribution of economic activity may 
shed light on a range of other systems that have been investigated separately using complex networks. 

PACS numbers: 61.43.Hv, 89.75.Da, 89.75.Hc, 89.65.Gh 



Standard models of urban economics have difficul- 
ties describing economic growth and as a consequence 
they cannot reproduce important empirical observations, 
maybe most notably Zipf's Law[l|,|2j,|3|]. One of the ba- 
sic tenets of market economics is that economic growth 
arises as a result of mutually beneficial trade between 
specialized agents. We use this mechanism as the basis 
of an abstract network representation of economic activ- 
ity and trade on a surface. 

Casting the urban economy as a complex network puts 
a range of new tools at our disposal for analyzing this im- 
portant system 0, 0, |(J. The benefit goes beyond direct 
observations about land uses and land values since many 
other growing spatial systems of human origin grow as a 
direct response to demand from within the economy. Be- 
cause of this, a model that simulates the geographic dis- 
tribution of economic activity can also shed light on the 
workings of many other social, technical and economical 
systems - many of which has been studies as complex net- 
works in their own right such as the Internet, phone call 
graphs and human sexual contacts 0, IE HI & CZJ E3 • 
It is quite conceivable that reported power law distribu- 
tions in these and other such systems may stem more or 
less directly from near-optimal distribution of services to 
meet demand within the geographic distribution of the 
urban system. 

We take fixed-size non-overlapping land areas to be 
network nodes, and the trade relations between them 
to be connections, in the network model. The activity 
within a node is defined as the number of connections 
connecting it to the rest of the network. We hereby ob- 
tain an abstract representation that can capture hetero- 
geneity in producer specialization and spatiality at the 
same time. This allows us to introduce a basic economic 
unit, the Economically Correlated Pair (ECP), as an ab- 
stract expected future stream of trade between two ac- 
tivities. Per definition, an ECP has unit economic value 
to both its end-points - a reasonable simplification given 
the very large number of trade relations in an economy. 



The simulations we use for the results in this Letter 
start from a lattice of undeveloped land areas and pro- 
ceeds to simulate the allocation of new urban activities 
and their trade relations. We compare statistics of simu- 
lated node degrees with empirically observed land values. 
The comparison assumes a linear relationship between 
the amount of money generated and the value of land. 
This relationship follows from our definition of ECPs and 
from i) market pricing of goods and services and ii) the 
connection between trade benefits and land value. 

i) Market pricing of commodities provides an adap- 
tive measure that allows us to compare the activities 
that generate them. Hence, on average, an edge con- 
tributes identically to the value of both nodes to which 
it connects. This contribution is exactly our definition 
of activity, which implies that the degree of a node is 
proportional to its benefits due to trade. 

ii) This connection consists of two proportionalities. 
For a node i we have 
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where Vi is the value of the corresponding land area, rj 
is the bid-rent 0| i an d %i is the total trade benefits as 
outlined above. Capitalizing periodic rent income from 
the site % g ives land value Vi = Q where i is the inter- 
est rate[l3|- The second proportionality is a weak form 
of the leftover principle from urban economics, which 
states that, in a competitive land market, rent equals 
the amount of money left after all expenses (except rent) 
are paid. This amount of money equals the sum of all 
trade benefits at the site. For our results it is sufficient 
that, on average, a certain proportion of each new unit 
of trade benefit goes to the landowner. 

We formulate the network model by dividing a ge- 
ographic area into N non-overlapping fixed-size cells 
{1,2,..., N} and taking these to be network nodes. The 
undirected connections between nodes are trade streams 
(ECPs) between activities in the cells, these can also con- 
nect a cell to itself and two nodes can be connected by 
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any number of mutual ECPs. The amount of profit (be- 
fore rent), generated in a node i, is per the definition of 
ECPs equal to mXj, where m is the value of an ECP and 
Xi is the degree. 

Network evolution is the process by which a network 
changes over time. New activity occurs in the form of ad- 
ditions of ECPs to the system. The growth mechanisms 
include balanced node addition and preferential attach- 
ment which gives a strong similarity to the Barabasi- 
Albert model, and a scale-free node distribution can be 
expected 

The network is initialized by the pairwise connection 
of no (even) spatially uniformly distributed nodes. Addi- 
tion of ECPs are separated into primary and secondary 
growth, corresponding to the placing of first and second 
edge end-points in the network. The placement of an 
end-point is done either preferentially (per-activity) or 
uniformly (per- node) ; corresponding to in situ expansion 
or establishment of new activity, respectively. Expansion 
of node activity can correspond to existing activities be- 
coming more space efficient as well as them being out-bid 
and replaced by more efficient activities. 

Preferential growth corresponds to trade with an ex- 
isting activity. Activities are considered to be average 
activities and thus they all interact to the same extent 
with each other. If we assume constant fractions q\ and 
<72 for uniform and preferential growth respectively, with 
Qi + 92 = 1, the probability of selecting a node i prefer- 
entially in primary growth is 
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Uniform growth corresponds to average activities select- 
ing sites in a competitive land market. The probability 
of selecting a node % uniformly as a primary effect is 
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where 
and S 
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is the number of developed nodes at time t 
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erwise. The meanings of n 

analogous to and 5^ D \ with P referring to perime- 
ter nodes and E referring to external nodes. Perime- 
ter nodes are all undeveloped nodes that are adjacent 
to a developed node and external nodes are undeveloped 
nodes without any developed neighbors. The parameter 
b controls the fraction of perimeter nodes that are open 
for development and e controls how many external nodes 
that are open for development compared to the number 
of perimeter nodes. The interpretation of this is that ex- 
ternal growth is identical to perimeter growth but takes 
place on ambient infrastructure not explicitly represented 
in the model. Infrastructure tends to grow as the urban 



system grows which makes it reasonable to assume that 
e is constant over time. 

The probability of secondary preferential growth at site 
i as a consequence of primary growth at site j is 
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and for secondary uniform growth, it is 
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Dij is a matrix of site-to-site interaction strengths, which 
decay with transportation costs. We have used Dij — 
(1 + cd(i,j))~ a , where d(i,j) is the Euclidean distance 
between sites i and j. The non- negative parameters c 
and a controls the impact of spatiality. 

To simplify the analysis of the model it is useful to 
assume that development of new land (addition of an ac- 
tive node) takes place at a constant rate qA (compared 
to other types of growth, not as a function of physical 
time). To motivate the assumption, let us consider the 
growth of developed clusters and their perimeters. It is 
true in simulations of the model, and it can be empiri- 
cally verified (see Figures^; and^l), that the cluster area 
distribution is close to a simple power-law with density 
function f(A) ~ , and that the relation between clus- 
ter perimeter P and cluster area A has the form P ~ A x , 
with A < 1. From this we observe that for the entire sys- 
tem of clusters wc have 
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assuming that > 2 and A < 1. We define q[ as the 
fraction of activity increments that occur uniformly on 
developed nodes, and Eq. |J5J gives 
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which, because of Eq. @, is approximately constant. 
This means that the rate of node activation, qA = 1 — 
q'i — 92, also can be considered constant. 

If we assume that the spatial bias between sites is 
small [li). the time evolution of expected activity on a 
developed site i follows the equation 



I x(t) 
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which is solved by the continuous-time method intro- 
duced by Barabasi et al . After sufficiently long time 
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FIG. 1: Diagrams (a), (b), and (c) show double-logarithmic histograms with exponentially binned empirical (x) and simulated 
(□) observables: (a) - land value per cell, (b) - aggregated cluster land value, (c) - cluster area. The lower part of the empirical 
distribution in (a) corresponds to land values below our urban land- use threshold m = 75kSEK/Ha, and is not described by 
our model. The exponents of the power law distributions in (a), (b) and (c) are roughly 2.1, 1.8 and 2.3 respective. In (d) is 
shown empirical (broad boxes) and simulated (thin boxes) results for cluster area plotted against exponentially binned cluster 
perimeters. The vertical interval of the boxes contains 90% of the perimeters in the bin. The reference line has a slope of 0.7. 
Part of the simulated configuration used for the statistics shown in the other figures is shown in (e). Diagram (f) shows cluster 
populations plotted against exponentially binned cluster prices. The vertical interval of the boxes contain 90% of the cluster 
prices in the corresponding bin and the crosses indicate the medians of the prices in the bins. The reference line has a slope of 
1, which indicates a linear relationship. The comparison was carried out by first identifying geographical clusters of high land 
values and then comparing their accumulated population content with their accumulated value. 



the degree distribution approaches the form P[xi = x] ~ 
(a; + B)-"t, with 



92(1 - q[ - 92) 

and 

7=1+—. (10) 

92 

Most of the parameters in the model can be readily 
estimated from empirical data. If n t is the number 
of clusters, we have e = n[ C ^ /n\ D \ The exponent 7 
in the distribution of activity gives information about 
the size of q^ (and q\) via Eq. (|10f) . The perimeter 

(D) 

parameter b can be determined by qA = *x-(t) an< ^ 

( u^\~ X 
b = qA I (31 -gA)(l + e)^fsy J • 

In principle, these parameters are not empirical in na- 
ture, they are aggregated measures of fundamental prop- 



erties of the current economic system (and region) of in- 
terest. The spatial parameters c and a reflect the statis- 
tics of transport characteristics of the economic configu- 
ration. They are not as easily estimated from data as the 
other parameters but results seem robust to their exact 
values and the functional form of Djj . 

The data used for all empirical results is based on a 
database delivered by Sweden Statistics that covers es- 
timations of the market value of all land in Sweden (2.9 
million data points). Cluster measurements was obtained 
by identifying clusters using land value as the threshold 
parameter (75 kSEK/Ha) and then comparing accumu- 
lated price and population within the obtained areas. To 
identify clusters in simulated and empirical data we use 
a computer program that masks away all data points 
below a threshold value and then treats all contiguous 
(8-cell neighborhood) areas as clusters. 

Once the frequency distribution of land value per unit 
area is correct, the challenge in realizing the cluster level 
is to faithfully capture the spatial sorting of the node 
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degrees. Using the model presented in this Letter we 
obtain excellent agreement between simulated and em- 
pirical statistics for a range of non-trivial higher-order 
structures in successive orders of upward causation (cells 
to clusters), see Figures JTJ)) through (TJa). These results 
are not trivial consequences of the land price distribution; 
it is perfectly possible to arrange the developed cells into 
any system of clusters. The same is true for the relation 
between cluster area and perimeter. Also, in Figure (JJJ) 
we demonstrate that geographical distribution of pop- 
ulation and land value can be observed interchangeably. 
This is a likely example of how an earlier studied spatially 
growing system of human origin likely is just a response 
to the growth of the economic network. The growth of 
the latter is defined by activity inter-dependency and 
transportation costs - a more fundamental candidate for 
a generating process than more vag ue postulations about 
why and when people re-locate [l?]]. 

All figures reflect data from the same Monte Carlo sim- 
ulation run. Figure JU) shows a part of the spatial con- 
figuration from this run. The parameters for the reported 
simulation, m = 75kSEK (10SEK« 1USD), b = 0.15 and 
e = 0.25 have been estimated from empirical data. By 
studying the exponent in the empirical distribution of 
land values, an initial estimate of q% =0.1 was obtained. 
For the reported results qi was adjusted to 0.2. The spa- 
tial parameters were c = 0.1 and a = 2. Investigation of 
sensitivity shows that exponents and proportions change 
slowly and smoothly with all parameters. A square grid 
of 2000 x 2000 cells was used and the number of iterations 
was 10 6 . 

Macroscopic models, such as Simon's Model and later 
derivations, are based on the rather unrealistic assump- 
tion that individual cities would react in unison to exter- 
nal events [l^, Lit l2i|- Although this may in some sense 
be trueQ, such models hold limited explanatory and pre- 
dictive power since they can not be mapped to a realistic 
situation where the internal structure of cities is impor- 
tant or when clusters interfere and coagulate. By making 
a similar assumption for small fixes-size land areas our 
representation remains valid regardless of the dynamics 
of higher-order structures such as cities. 

The relation between pattern (probability distribu- 
tions etc.) and process (microdynamics) is not one-to- 
one. Rather, for each observed macroscopic pattern there 
is likely a whole class of processes capable of generating 
it. Apart from arguing for the network model ontology 
on the micro level there is a number of implications that 
remain to test against empirical data. Among these are 
the time evolution of land values, explicit measurement 
of the trade network and making studies similar to the 
one presented here for other countries. 
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